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F. : Abstract 



A light-front treatment for finite nuclei is developed from a relativistic 
effective Lagrangian (QHD1) involving nucleons, scalar mesons and vector 
mesons. We show that the necessary variational principle is a constrained 
one which fixes the expectation value of the total momentum operator P + 
to be the same as that for P . This is the same as minimizing the sum 
of the total momentum operators: P~ + P + . We obtain a new light-front 
version of the equation that defines the single nucleon modes. The solutions 
of this equation are approximately a non-trivial phase factor times certain 
solutions of the usual equal-time Dirac equation. The ground state wave 
function is treated as a meson-nucleon Fock state, and the meson fields are 
treated as expectation values of field operators in that ground state. The 
resulting equations for these expectation values are shown to be closely related 
to the usual meson field equations. A new numerical technique to solve the 
self-consistent field equations is introduced and applied to 16 and 40 Ca. The 



computed binding energies are essentially the same as for the usual equal- 
time theory. The nucleon plus momentum distribution (probability for a 
nucleon to have a given value of p + ) is obtained, and peaks for values of 
p + about seventy percent of the nucleon mass. The mesonic component of 
the ground state wave function is used to determine the scalar and vector 
meson momentum distribution functions, with a result that the vector mesons 
carry about thirty percent of the nuclear plus-momentum. The vector meson 
momentum distribution becomes more concentrated at p + = as A increases. 

I. INTRODUCTION 

The purpose of this paper is to derive the light-front formalism necessary to compute the 
properties of finite nuclei. Nuclear properties are very well handled within the existing con- 
ventional nuclear theory, so it behooves us to explain why we are embarking on this project. 
Our motivation is that understanding experiments involving high energy nuclear reactions 
seems to require that light-front dynamics and light cone variables be used. Consider the 
EMC experiment [|J, which showed that there is a significant difference between the parton 
distributions of free nucleons and nucleons in a nucleus. This difference can interpreted as a 
shift in the momentum distribution of valence quarks towards smaller values of the Bjorken 
variable xbj- The Bjorken variable is a ratio of the plus-momentum k + = k° + k 3 of a quark 
to that of the target. Thus light cone variables are relevant. If one uses k + as a momentum 
variable the corresponding canonical spatial variable is x~ = x° — x 3 and the time variable 
is x + = x° + x 3 . 

It is important to realize that the use of light-front dynamics is not limited to quarks 
within the nucleon - - it also applies to nucleons within the nucleus. This formalism is 
useful whenever the momentum of initial or final state nucleons is large compared to their 
mass ||. In particular, it can be used for (e, e'p) and (p,2p) reactions. If one uses light- 
front variables for nucleons in a nucleus, it is also necessary to maintain consistency with 
the information derived previously using conventional nuclear dynamics. This provides the 
technical challenge which we address in the present manuscript. 

The conventional equal-time approach to nuclear structure physics provides an excellent 
framework, so it is worthwhile to introduce the light-front variables and describe the expected 
advantages in a general way. The use of the light-cone variables can be obtained using a 
simple argument based on kinematics 0. Suppose the virtual photon that is absorbed by a 
fermion at a space-time point (z±, t\). The fermion then starts to move at high momentum 
and nearly the speed of light and emits the photon at another space time point (z2, £2)- 
In between the two times, the wave function of the entire system has undergone a time 
evolution given by the complicated operator e~ %H< ^ t2 ~ tl \ But we have z\ + ct\ = zi + cti-, 
if the z— axis is opposite to the direction of the virtual photon. The two scattering events 
occur at different times, but at the same value of x + = z + ct. Thus if we use x + as a time 
variable, no time evolution factor appears. The net result is that the cross section involves 
light-like correlation functions which involve field operators evaluated at the same light-front 
time: x + = (see for example the reviews RH]). Thus it is a specific and general feature of 
the light-front wave approach that knowing only the ground state wave function is sufficient 



for computing the distribution functions. 

Let us review the salient features of the basic idea that using the light-front approach 
leads to a simplified treatment. To be specific, consider high energy electron scattering 
from nucleons in nuclei. The key ingredient in the light-front simplification is to realize 
the main difference between the two formalisms. In the equal-time formalism, sums over 
intermediate states are taken over eigenstates of the Hamiltonian, P°. The usual three 
momentum is conserved, but energy is not conserved in intermediate states. In the light- 
front approach one sums over eigenstates of the minus component of the total momentum 
operator. The value of P~ is not conserved in intermediate state sums, and the values of P + 
and P 1 - are conserved. This is especially convenient for high energy reactions, in which the 
plus-component is the largest component of the momentum for each projectile or ejectile. 

The advantage of using P~ as an "energy" variable can be easily described. Let the 
four-momentum q of the exchanged virtual photon be given by (u, 0,0, — \JQ 2 + z/ 2 ), with 
Q 2 = —q 2 , and Q 2 and v 2 are both very large but Q 2 jv is finite (the Bjorken limit). In this 
case it is worthwhile to use the light-cone variables q ± = q°±q 3 in which q + ps Q 2 j2v = Mx, 
q~ « 2v — Q 2 j2v, so that q~ ^> q + . Here M is the mass of a nucleon and x is the Bjorken 
variable. We shall neglect q~ in comparison to q + , noting that corrections to this can be 
handled in a systematic fashion. Then the schematic form of the scattering cross section for 
e + A — > e' + (A — 1)/ + p, where / represents the final nuclear eigenstate of P~ , and p the 
four-momentum of the final proton, is given by 

da ^^ J ^ J d 4 pS(p 2 - M^^iq + p.-pj -p)\(pj \ J(q) \t) \ 2 . (1.1) 

j ■' 

Here the operator J(q) is a schematic representation of the electromagnetic current. Per- 
forming the four- dimensional integral over p leads to the expression 

da ~ £ / d ^I^J- b fa - pf + q f - m 2 ) | (p, / | J {q) | i) \ 2 . (1.2) 

f J Pf 

The argument of the delta function (p^ — pj + q) 2 — M 2 pa — Q 2 + 2q~(pi — Pf) + ■ Thus 
we see that p7 does not appear in the argument of the delta function, or anywhere else, 
so that we can replace the sum over intermediate states by unity. In the usual equal-time 
representation, one finds the argument of the delta function to be — Q 2 + 2v(Ei — Ef). The 
energy of the final state appears, and one can not do the sum. 
To proceed further in this schematic approach we take 

J(q) = J d 3 kb{ +q b k , (1.3) 

where b is a nucleon destruction operator and V = (V + , V±). It is useful to define ps = 
Pi — pf because 

p + B = Q 2 /2v = Mx, (1.4) 

as demanded by the delta function. Then one can re-express Eq. (|1.2| ) as 



da ~ d 2 p B± (i | b j pB b PB \ i) = d 2 p B j_n(Mx,p BA _), (1.5) 

where n(Mx,pB±) is the probability for a nucleoli in the ground state to have a momentum 
(Mx,pb±)- Integration in Eq. (|1.5|) leads to 

a ~ f d 2 p ± n(Mx,p±) = f{Mx), (1.6) 

with f(Mx) as the probability for a nucleon in the ground state to have a plus momentum 
of Ma;. 

The quantity f(Mx) has been a widely used prescription UP for handling the light-front 
in a simple way. The variable Mx is replaced by M — e a + k 3 , in which the label a denotes 
a shell-model orbital <p a of binding energy e a . Then 

fup = Y, n * [d 2 p± [dp 3 \0 a (p 3 ,p ± )\ 2 5(M-e a +p 3 -Mx), (1.7) 

in which n a is an occupation probability. The validity of this prescription, which rests on 
a reasonable assumption, is rather suspect because the variable p + = Mx is a kinematic 
variable, unrelated to discrete eigenvalues of a wave equation. One of the main purposes of 
the present paper is to see if anything like this prescription emerges from our calculations. 
We shall see that Eq. ( |1.7|) is not obtained, if a vector potential is a significant part of the 
nuclear mean field. 

It is useful to discuss the relation with ^-scaling ||. The arguments that the cross 
section depends on a plus-momentum distribution are well known when used for quarks in 
a nucleon, but they also apply to nucleons in a nucleus |J. Ji and Filippone || showed that 
the y scaling function F(y) extracted in quasi-elastic electron scattering on nuclei is actually 
the light cone plus-momentum distribution function for nucleons in the nucleus. It is useful 
to use a relativistic form of the variable y J7| in which 

y = -q 3 + u + E s + M, (1.8) 

as both q 3 and v are large in magnitude, and E s is the single nucleon separation energy. But 

(q 3 ) 2 -u 2 
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1 + %^, so that 



Mx = M + E s -y = M A y A . (1.9) 

Here y A is a new y-scaling variable. This means that according to Eq. ( |1.4| ) 

p+ = M A y A /M w Ay A , (1.10) 

so that a measurement of a determines the probability that the struck nucleon has a plus- 
momentum of Ay A . This probability also enters in convolution model calculations of nuclear 
deep inelastic scattering. 

The use of light-front dynamics to compute nuclear wave functions should allow us to 
compute F(y) from first principles. Furthermore, we claim that using light-front dynamics 
incorporates the experimentally relevant kinematics from the beginning, and therefore is the 



most efficient way to compute the cross sections for nuclear deep inelastic scattering and 
nuclear quasi-elastic scattering. 

It is worthwhile to review some of the features of the EMC effect [[M] . The key experi- 
mental result is the suppression of the structure function for x ~ 0.5. This means that the 
valence quarks of bound nucleons carry less plus-momentum than those of free nucleons. 
One way to understand this result is to postulate that mesons carry a larger fraction of the 
plus-momentum in the nucleus than in free space. While such a model explains the shift in 
the valence distribution, one obtains at the same time a meson (i.e. anti-quark) distribution 
in the nucleus, which is strongly enhanced compared to free nucleons and which should be 
observable in Drell-Yan experiments ||. However, no such enhancement has been observed 
experimentally ||, and the implications are analyzed in Ref. [10 |. 



The use of light-front dynamics allows us to compute the necessary nuclear meson dis- 
tribution functions using variables which are experimentally relevant. The need for a com- 
putation of such functions in a manner consistent with generally known properties of nuclei 



led one of us to attempt to construct a light front treatment of nuclear physics [[iT|] . These 
calculations, using a Lagrangian in which Dirac nucleons are coupled to massive scalar and 
vector mesons |l^j , treated the example of infinite nuclear matter within the mean field 
approximation. In this case, the meson fields are constants in both space and time and the 
momentum distribution has support only at k + = 0. Such a distribution would not be acces- 
sible experimentally, so that the suppression of the plus-momentum of valence quarks would 
not imply the existence of a corresponding testable enhancement of anti-quarks. However, 
it is necessary to ask if the result is only a artifact of the infinite nuclear size and of the 
mean field approximation. The present paper is an attempt to handle finite-sized nuclei 
using light-front dynamics. 

A. Recovery of rotational invariance 

It is worthwhile to discuss, in a general way, how it is that we are able to find spectra 
which have the correct number of degenerate states. Let us imagine that we try to determine 
eigenstates of a LF Hamiltonian by means of a variational calculation. Simply minimizing 
the LF energy obviously leads to nonsensical results since the LF energy scales like the 
inverse of the LF momentum. Even if one has only a poor ansatz for the intrinsic wave 
function, one can easily reach zero energy by letting the overall momentum scale to infinity! 
However, this problem is avoided by performing a constrained variation, in which the total 
LF momentum is fixed by including a Lagrange multiplier term proportional to the total 
momentum in the LF Hamiltonian. Note that this is not a problem if one is able to use a 
Fock space basis in which the total plus and _L momentum of each component are fixed. In 
calculations involving many particles, the Fock state approach cannot be used in practical 
calculations — instead one uses a mean field in which each particle moves in an "external" 
potential. In this case the total momentum is not fixed, and a Lagrange multiplier term 
needs to be included in order to avoid solutions with infinite LF momentum. 

In order to fix this potential problem with "runaway solutions" (P + — » oo) to variational 
calculations for LF Hamiltonians, any term proportional to P + would suffice. However, 
by setting the coefficient for the term proportional to P + equal to one, i.e. minimizing 



P +P + , one automatically guarantees that P + = P (or P 3 = 0). The reason is that, using 
covariance, P~ has eigenvalues of the form P~ = n p+ x , i.e. it scales like 1/P + . Therefore, 
when one minimizes P~ +P + with respect to P + , the minimum occurs for P + = J M% + P]_, 

which yields P~ = n P + ± = yM 2 + P± as well. This "equipartition" between P + and P~ 
thus arises since the two operators scale in exactly opposite ways under longitudinal boosts. 
Note that this is quite analogous to the nonrelativistic harmonic oscillator where, under 
scale transformations, potential and kinetic energy scale in opposite ways, resulting in the 
equipartition between potential and kinetic energy. 

The net result is that we minimize the sum of P + + P~ . The need to include the plus- 
momentum can also be seen in a simple example. Consider a nucleus of A nucleons of 
momentum Pj[ = Ma, Pal = 0, which consists of a nucleon of momentum {p + ,p±), and a 
residual (A — 1) nucleon system which must have momentum (P4 — p + , —p±). The kinetic 
energy K is given by the expression 

p1± M 1 A + Ml_ (1U) 



rvj 



In the second expression, one is tempted to neglect the term p + in comparison with P\ 
Ma- This would be a mistake. Instead make the expansion 

pj + M 2 M\_ x ( p+\ 

p + n \ n) 

v 2 + M 2 

~^- +p + + M A -i, (1.12) 

p+ 

because for large A, M\_ x jP\ f» 1. For free particles, of ordinary three momentum p one 
has E 2 (p) = p 2 + m 2 and p + = E(p) + p 3 , so that 

K « ( E2 (P)-(^) 2 ) + E {p) + p 3 + M A _! = 2E(p) + M A -l (1.13) 

We see that i^ depends only on the magnitude of a three-momentum and rotational invari- 
ance is restored. The physical mechanism of this restoration is the inclusion of the recoil 
kinetic energy of the residual nucleus. 



B. Outline 

The organization of the paper is as follows. The light-front quantization for our chosen 
Lagrangian is is presented in Sec. II. This quantization is applied, along with a constrained 
minimization of the expectation value of P , to derive a light-front version of mean field 
theory in Sec. III. We obtain a new light version of the equation that defines the single 
nucleon modes. The solutions of this equation are approximately a non-trivial phase factor 
times the solutions of the usual equal-time ET Dirac equation. The consequences of this 
phase factor are discussed. 



The meson fields are treated as expectation values of operators. The equations for these 
expectation values are closely related to the meson field equations appearing in the usual 
treatment of the Walecka model. However, the mesonic Fock space is accessible in our 
formalism. Our nucleon mode equation is simplified by the use of a two-component spinor 
formalism ||13| , and by an angular momentum reduction in Sec. IV. The numerical aspects 
are discussed in App. [A]. The binding energies, nucleon and meson distributions for 16 and 
40 Ca are presented in Sec. V. A concluding discussion appears in Sec. VI. Numerical details 
of how we evaluate the momentum distributions are given in App. ||. A brief discussion of 
some of the results can be found in Ref. [H|. A related set of solutions of some toy model 



problems and a heuristic derivation of our nucleon mode equation will appear in a separate 
paper |15 . 



II. LIGHT-FRONT QUANTIZATION 

We start with a model in which the nuclear constituents are nucleons ip (or ?/>'), scalar 
mesons <fi and vector mesons V^. The Lagrangian £ is given by 

C = \{d^3^ - m 2 s 4> 2 ) - \v^V, v + \m\y^ + fi (V (| d, -g v V„)-M- g,<?) ij/, 

(2.1) 

where the bare masses of the nucleon, scalar and vector mesons are given by M,m s , m v , 
and V^ u = d^V u — d u V^. We ignore pions here. 
The field equations are given by: 

j -(id- g v V)if/ = (M + grfW, (2.2) 

d,V^ + mlV v = s^W (2-3) 

9„^ + m^ = -^y. (2.4) 



The next step is obtain the light-front Hamiltonian (P~) JTBJ as a sum of a free, non- 
interacting and a set of terms containing all of the interactions. This is accomplished by 
separating the independent and dependent degrees of freedom in the usual way fTTPf and 
then using the energy momentum tensor. Consider the nucleons. Although described by 
four- component spinors, these fields have only two independent degrees of freedom. The 
light-front formalism allows a convenient separation of dependent and independent variables 
via the projection operators A± = |7°7 ± J18|jl9|| , with tp± = A±i/)'. The independent fermion 



degree of freedom is chosen to be ip' + , and one finds 

(id~ - g v V-)iP' + = (ol • (p ± - g v V ± ) + (3{M + g s <f>))<ip'_ 

{id + - g v V + )iP'_ = (ol • ( P± - g v V ± ) + (3(M + g s (f>W + . (2.5) 

The relation between ip'_ and if)' + is very complicated unless one may set the plus component 
of the vector field to zero fl?||. This is a matter of a choice of gauge for QED and QCD, but 



the non-zero mass of the vector meson prevents such a choice here. Instead, one simplifies 
the equation for if)'_ by |TE| , |2"0| | transforming the fermion field according to 



if/ = e~ i9vA{x) ip, d + A = V + . (2.6) 

This transformation leads to the replacement of Eq. ( |2.5| ) by 

(id- - g v V-)ip + = (ol ■ ( P± - g v V ± ) + (3(M + g s (j>))ip- 

td + ^ = (ai ■ (p± - g v V ± ) + /3(M + g s <f>))ip + , (2.7) 

where 

d+V^ = d + V^ - d^V + . (2.8) 

Note that while it is V^ that enters in the nucleon field equations, it is V^ that enters in 
the meson field equations. 

The scalar field can be expressed in terms of creation and destruction operators: 



d 2 k ± dk+6(k H 

IX) 



(27r) 3 / 2 v / 2AF 
where 



a(fc)e- Jfe - x + a T (fc)e , ( 2 -9) 



k ■ x = -{k + x ) — k± ■ x± = k ■ x, (2-10) 

and the fields and their derivatives with respect to x + are evaluated at x + = 0. This notation 
is used through out this work. The consequence is that the energy momentum tensor T^ v 
does not depend on x + . In the above expansion (and in the expansions for any of our fields) 

the particles are on the mass-shell. Here k~ = \+ s . The theta function restricts k + to 
positive values. The commutation relations are 

[o(fc), a f (fc')] = S(k± - k' ± )5{k + - k' + ), (2.11) 

with [a(k),a(k')] = 0. It is useful to define 

<5 (2 ' +) (fc - k') = 5(k ± - k' ± )5(k + - k' + ). (2.12) 

The expression for the vector meson field operator is 

V "W = f fo^r^ £ ^(k, w) [a(fc, W )e"*- + fl t(fc, W )e*-] , (2.13) 

J (27r) 3 / 2 V2k+ w=1)3 L J 

where the polarization vectors are the usual ones: 

fc%(fc, u) = 0, e M (fc, w)e^(fc, u/) = -<W, 

£ ^(fc )W )e"(* )W )=-(^-^ T ). (2.14) 

w=l,3 m v 

Once again the four momenta are on-shell, with k~ = ± k+ - . The commutation relations 
are 
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[a(k,u),a\k',u')] =6^ 2 '+\k-k'), (2.15) 

with [a(k,cj), a(k',uj')] = 0, and lead to commutation relations amongst the field operators 



that are the same as in Ref. 20 



We also need the eigenmode expansion for V^. This is given by 



V(x) 



d 2 k ± dk + 6(k- 



E **(*■ 



w 



[k,u;)e- ik - x + aHk,Lu)e 



ik-x 



where, using Eqs.(|2~8|) and ( |2.13| ), the polarization vectors e M (fc,u;) are 

k^ 
e^(k,uo) = e"(Jfe,w) - —e + (k,u). 

k + 

Note that 

k v k^ 

^ e»(k,u)e"(k,U,) = -tf>>-g+l>—-g+> — ). 

w=l,3 

Then we may construct the total four-momentum operator from 
p^ = - fdx-d 2 x ± T + v{x + = 0,x~,x x ), 
with (as usual) 



(2.16) 



(2.17) 



(2.18) 



T" u = -g^C + J2 



9C SVM 



(2.19) 



(2.20) 



in which the degrees of freedom are labelled by <p r . We need T ++ and T + , which are 

T++ = d + (j)d + (j) + V ik V ik + m 2 v V + V + + 2iplid + ip + , (2.21) 

and 



T+- = V ±( j) • V ± + m\<\) 2 + ^(V + -) 2 + ]-V kl V kl + m^V fe 



+V h± ■ (P± ~ 9vV~) +M + g 8 <f>) i,. 



(2.22) 



This form is still not useful for calculations because the constrained field ip- contains 
interactions. We follow Refs. [H|21] in expressing i\)_ as a sum of terms, one £_ whose 
relation with ijj + is free of interactions, the other r/_ containing the interactions. That is, 
rewrite the second of Eq. 



v- 



as 
1 

1 
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[a± ■ p± + /3M)il>+ 
{-oc± ■ g v V± + (3g s <p)4> + . 



(2.23) 



Furthermore, define £+(x) = ip+(x), so that 



ip{x) = £(x) + ?7_(x) 



(2.24) 



where £(x) = £_(x) + £+(x). This separates the dependent and independent parts of ip- 

One needs to make a similar treatment for the vector meson fields. The operator V + ~ , 
is determined by 



V 



-+ 



fr< 



g v J + - m 2 v V + - diV 1 



(2.25) 



Part of this operator is determined by a constraint equation, because the independent vari- 
ables are V + and V l+ . To see this examine Eq ( |2.25| ), and make a definition 



V 



where 



CO 



oj 



d+ J 



(2.26) 



^2.27) 



The sum of the last term of Eq' ( p. 22 ) and the terms involving u + ~ is the interaction density. 
Then one may use Eqs. ( |2.22|) , (|2.24|) , and ( |2.26| ) to rewrite the P~ as a sum of different 
terms, with 



Pqn = § / d 2 x ± dx- e (7.L ■ P± + M) f , 



and the interactions 



with 



Pj = Vi +v 2 +v 3 , 



Vl 



j d 2 x ± dx £ (g v i ■ V + g s ^j £, 



(2.28) 



(2.29) 



^2.30) 



and 



V2 



(2.31) 



J d 2 x ± dx £ (-g v l ■ V + g s 4>) —^ {~9vl • V + g a <f>) f , 

V3 = ~i J d2x±dx ~ J dy i e ( x ~ ~ y ^ £+(^> ^±)7 + c+(i/r, «±) 

x y dyj c(x~ - Z/2 _ )d(^> a;±)7 + ^+(^> ®±), (2-32) 

where e(x) = #(:r) — 6(—x). The term t>i accounts the emission or absorption of a single 
vector or scalar meson. The term t> 2 includes contact terms in which there is propagation 
of an instantaneous fermion. The term v 3 accounts for the propagation of an instantaneous 
vector meson. 

Our variational procedure will involve the independent fields ip + , so we need to express 
the interactions Pq N , t?i j2 in terms of £ + . A bit of Dirac algebra shows that 
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P N = PqN + V!+V 2 



d 2 xJ^-&[2g v V- 



+ (a ± ■ (p± - g v V ± ) + {3{M + g,<j>)) — (a± ■ (p± - g v V x ) + /3(M + g,<f>j)]t+. (2.33) 

It is worthwhile to define the contributions to P arising from the mesonic terms as P^ 
and P^. Then one may use Eqs. fl2.22|) and ( p.21|) along with the field expansions to obtain 



p- = - d 2 x±dx~ (V ± ■ V ± + m 2 s (j) 2 

k 2 , +m 2 



d 2 k ± dk + e(k + )a\k)a(k)'- 



k+ 



(2.34) 



P+ = f d 2 k ± dk + e{k + )a j {k)a{k)k H 



(2.35) 



and 



P~ = J2 I d 2 k ± dk + 6{k A 



2 i ™2 



k i + m: 



w=l,3 



k+ 



-a Jf (k,u)a(k,uj) +V3, 



P+ = Y, I d 2 k ± dk + e{k + )k + a\k,uj)a{k,u}). 

. . — 1 Q J 



w=l,3 



(2.36) 



(2.37) 



The term V3 is the vector-meson instantaneous term, and we include it together with the 
purely mesonic contribution to P~ because it is cancelled by part of that contribution. 
Thus, our result for the total minus-momentum operator is 



and for the plus-momentum 



P — P N + P s + P v , 



P+ = P+ + P+ + P V + , 



where from Eq. ( |2.21 



P+ = Jd 2 x ± ^2&d+i 



(2.38) 



(2.39) 



(2.40) 



III. MEAN FIELD THEORY 



The light-front Schroedinger equation for the complete nuclear ground-state wave func- 
tion I \l/) is 



P-\V) = M A \ *). 



(3.1) 
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We choose to work in the nuclear rest frame so that we also need 

p+ | #) = m a | *). (3.2) 

We want to use a variational principle. One might think that one may simply minimize the 
expectation value of P~ , but this makes no sense because P + P~ = M\ when acting on the 
wave function. One would get a zero of P~ for an infinite value of P + . As explained in the 
Introduction, one must minimize the expectation value of P~ subject to the condition that 
the expectation value of P + is equal to the expectation value of P~ . This is the same as 
minimizing the average of P~ and P + , which is the rest-frame energy of the entire system. 
To this end we define a light-front Hamiltonian 

H LF = 1 -(P+ + P-). (3.3) 

We stress that H FF is not usual the Hamiltonian, because the light-front quantization is 
used to define all of the operators that enter. 

The wave function | \I/) consists of a Slater determinant of nucleon fields | $) times a 
mesonic portion 

| #) =| $)<8> | mesons), (3.4) 

and the mean field approximation is characterized by the replacements 

0^ (* | | *) 
V* -»• <* | V" I *). (3.5) 

We shall derive the meson field equations, and then determine the nucleon modes using a 
variational principle. 

A. Meson field equations 

We shall go through the derivation of the equation for the expectation value of <f>{x) in 
a detailed fashion. Consider the quantity H LF a{k) | \1>), and use commutators to obtain 

H LF a(k) | *> = [H LF , a(k)} | *> + M A a(k) | *). (3.6) 



The operators P s of Eqs. ( |2.34j ) and ( |2.35| ) and the standard commutation relations allow 
one to obtain 

[H LF , a(k)} = - kl + k ~[l + m2s a(k) + ^L (3.7) 

1 V n 2k+ K J (2vr) 3 / 2 v / 2fc+ 

where 3/ f 2 4= is the commutator of the interaction, Eq. ( |2.29| ), between the scalar meson 
and the nucleon: 

J{k) 1 



(27r) 3 /2 v ^T 2 
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[P 7 -,a((fc)]. (3.8) 



We use Eqs. ( [2.30|) -( p.32|) , and take the commutator of the interactions Vi with a(k). 



Then re-express the results in terms of £ to obtain 

J(k) = -\g s j d 2 x ± dx~ e ik ^(x)^{x). (3.9) 

Take the overlap of Eq. (|3.6| ) with (\1> | to find 



V ' V ; ' 7 2fc+ (2tt) 3 / 2 v / 2T+ V ; 

Multiply the above Eq. ( p.lOj ) by a factor /2^\z/i ^~' lk ' x ■ Then add the result of that operation 
to its complex conjugate. The integral of the resulting equation over all k± and positive 
values of k + and using the field expansion ( p.9|) leads to the result 

(<p | J <? k ±dk+9(k + ) (J(k)e~ ih - X + f(k)e +ih ' x ) | ¥). (3.11) 

The evaluation of the right-hand-side of Eq. ( p.ll|) proceeds by using Eq. ( |3~9D and its 
complex conjugate. The combination of those two terms allows one to remove the factor 
9{k + ) and obtain a delta function from the momentum integral. That |fc + appears in the 
exponential leads to the removal of the factor ~ of Eq. ( |3.9| ). One can also change variables 
using 

x = (z,x±). (3.12) 



2 

The minus sign enters to remove the minus sign between the two terms of the factor k ■ x 
in Eq. (|2.10 ). Then one may re-define the operator — V^ — [2-^-\ appearing in Eq. fl3.11|) 



as —V 2 . Note that we previously J22| obtained the above relation |3.12j simply by examining 



the space-time diagram for a static source (independent of x°). The net result is that 

(-V 2 + m 2 ) (* | <j)(x) | *) = -g B (V | 4>(x)iP(x) \ *), (3.13) 

which has the same form as the equation in the usual equal-time formulation. Note that the 
right hand side of Eq. ( |3.13| ) should be a function of \x\ for the spherical nuclei of our present 
concern. Our formalism for the nucleon fields uses x± and x~ as independent variables, so 
that obtaining numerically scalar and vector nucleon densities that depend only x\ + (x~ /2) 2 
will provide a central, vital test of our procedures and mean field theory. Assuming for the 
moment that this occurs, the scalar field (\I/ | (j){x) \ \&) will depend only |cc| according to 



(ra 



We stress that the use of Eq. (|3.12|) is merely a convenient way to simplify the calculation 



- using it allows us to treat the _L and minus spatial variables on the same footing, and 
to maintain explicit rotational invariance. We will obtain the mesonic plus-momentum 
distributions from the ground state expectation value of different operators. 
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The procedure of Eqs. ( |3.6| ) to ( |3.13|) can also be applied to the vector fields. The 
appearance of the barred vector potential makes it necessary to display certain steps. The 
starting point is to consider the expression HLpa(k) | \I/) and the interaction 

Using equations (|2.30|) - (|2.32|) , taking the commutator of the interactions Vi with a(k), and 



re-expressing the results in terms of £, leads to 

J(k, to) = --g v f d 2 x ± dx~ e ik - x £(x)-f ■ e(k, uj)£{x). (3.15) 

This, along with the other terms in the expression for ifu?a(A;, u) \ ^f), allows us to obtain 

at, I (u M^ k l + k +2 + ml (* I -lg v S<Px ± dx-e ik -*Z(x)'ye(k,u)Z(x) | g) 

(W a(K,oj) M/> — = ; — ^= . (3.1o) 

N ' V ' ' ' 2fc+ (27T) 3 /2 v ^fc+ V ' 

The field equation for the vector mesons V^ 1 is obtained by multiplying the above by e^ l (k, u), 
summing over u and performing standard manipulations. We need to know the quantity 

X"(k) = E7/(fc,w)6"(fc, W ) (3.17) 

The use of Eqs. ( gig) and ( fTTTD leads to 



X"(k) = -y + 28(i,, -)^ + ^7 + - (3.18) 

One makes familiar manipulations to obtain the result 

(-V 2 + m 2 v ) <* | V^x) | *> = g v {* | i{x)^i{x) | tt> + A", (3.19) 

with 

A"(«) = -(* I / |^ e*<— \>') ^ 7 + + 2*0*, -)^ £(*') I *•) (3.20) 

Note that (as in the derivation given above for 0) the variable fc + (confined to positive 
values) is replaced by the inclusion of the complex conjugate term (a T (fe, uj) by a variable fc 3 
which ranges from — oo to oo. We proceed by first assuming that 

(-V 2 + m 2 v ) (V | V(x) | tf) = g v (V | f(»)ye(aj) I *> = J", (3-21) 

which is to be verified by showing that Eqs. (|3.19| ) and (|3.21|) are consistent with the defining 



relation (|2.8|) Taking the difference between Eqs. ( [3.19D and ( p.21| ) and using the defining 
relation leads |2<| to a consistency requirement 

8»J+ = d + A^. (3.22) 

For fj, ^ — the above relation is verified by integration by parts in the expression (|3.20| ) for 
A M . If \x — — , one may use the definition Q2.8|) and that (\1/ | V + | \1/) does not depend on 
x + to see that 

dr{V\V- |*)=a-(*|V _ I*}. (3.23) 

Thus the validity of Eq. ( |3.21| ) is established. 
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B. Nucleon single-particle wave functions 

The mesonic field equations are given in the previous subsection. The equation for the 
nucleon modes are to be found using the procedure of minimizing P +P + with respect to the 
nucleon wave function, subject to the condition that the normalization of the independent 
fields remains fixed. The nucleon field operators enter only in the term P^ + Pjy, so that it 
is useful to define 

Hlf = \{Pn+ P$) ■ (3-24) 



The specific operator is obtained by using Eq. ( |2.33| ) to find 

dx 



/ft rp 
—d 2 x ± &H LF £ + , (3.25) 



where 

2H LF = id + + 2g v V~ 



«i • (p± - g v V±) + (3{M + g s <f))) — (a± ■ (p± - g v V ± ) + (5{M + g a <j>)) . (3.26) 



The potentials appearing in Eq. ( |3.26|) are independent of x + . This implies some sim- 



plifications: d+V~ = d + V- - d-V + = d+V~, so that V~ = V~, and (for i = 1,2) 



d+V % = d+V i - d % V + = -d i V + . Using the relation Q we find that V i = -9* A. 

The Slater determinant | $) is defined by allowing A nucleon states, denoted by the 
index a to be occupied. For our Slater determinant the constrained minimization is given 
by the equation 

sJd? x± *?L- (a | A + (n LF - P -A A+ | a) = 0, (3.27) 

where the quantities p~ are the Lagrange multiplication factors for each occupied orbital. 
The relation ( 3.27| ) leads immediately to our mode equation 



p-A + | a) = (id + + 2g v V~) A_ | a) 

+ (a ± ■ (p± - g v V ± ) + (3{M + g 8 <f>j) — (a ± ■ (p x - g v V ± ) + (3{M + g s <f>)) A + | a). (3.28) 

The operators a and (3 have non-zero values when appearing between A + and A_, but 
vanish when appearing between the two identical projection operators. Thus we may obtain 

A_ | a) as 

A- | a) = -^ (a ± ■ (p ± - g v V ± ) + (3{M + g 8 <j>)) A + | a), (3.29) 



or 



id+ 

i<9 + A_ | a) = (a± ■ (p± - g v V ± ) + j3{M + g 8 <f>)) A + | a). (3.30) 
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One may use Eq. (|3.30|) to rewrite Eq. ( |3.28|) as 
p-A + | a) = (id + + 2g v V~) A+ | a) 



a± ■ (p± ~ 9vV ± ) + (3{M + g s cj))) A_ | a) 



(3.31) 



Equations ( |3.31| ) and ( |3.30| ) are the essential results of this section. We have obtained the 
light-front version of the Hartree equations. 



C. Nuclear energy 

There are contributions to the expectation value of P + P + from the nucleons, scalar 
mesons, and vector mesons. The nucleonic term is given from the expectation value of the 
nucleonic part H LF ( |3.26| ). Taking the nuclear expectation value of H LF leads to a sum of 
matrix elements in the occupied states | a) . The use of the wave equation (|3.28|) leads to 
the result 



(* | H LF | *) = J2 



Vc 



(3.32) 



in which the sum is over a includes only occupied states. 

The contribution from the scalar mesons E s is given from the scalar meson terms of 
Eqs. (|2~34] ) and (g3|) by 



E M 



d 2 k±dk + 9{k^ 



k\ + m 2 s 
k+ 



+ k A 



(* | a\k)a(k) \ *). 



In our mean field approximation 

(* | a\k)a(k) | *) = |(# | a{k) \ *) | 2 , 



^3.33) 



(3.34) 



with the matrix element already known from Eq. ( p.lOj ). Then straightforward calculation 
leads to the result 



E x 



d 3 k 



2 J (2vr) 3 k 2 + ml 



(V | J(k) | *> 



(3.35) 



where J(k) is given by Eq. ( |3.9|) , the replacement ( p . 1 2|) is used, and as above k + is replaced 
by k 3 . The above expression is strikingly familiar --it is the result obtained in standard 
equal-time calculations. 

The vector meson contribution to the energy E v is defined as the one-half of the sum of 
the terms of P^ 1 of Eqs. ( |2.36|) and (|2.37|) . The calculation of P^ 1 is rather similar to the 
one just done for the scalar mesons. One uses the results (|3.15| ), (|3.16|) , and X M (|3.18| ). The 
effects of the instantaneous term v^ are cancelled by the non-7 M term of X M , so that we find 



E„ — — - 



d 3 k 



2 J (2vr) 3 k 2 + m 2 v 



J v (k) 



(3.36) 
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where 

J v (k) = g v (V | fd 3 xe ik - x ^(x)^(x) | *}. (3.37) 

The nuclear mass Ma is then given by 

occ - 

M A = Y,^ + E S + E V , (3.38) 

with expressions for each of the contributions given above. 

D. Relation with the equal-time formulation 

Our main results obtained using the mean field approximation and including the recoil 
of the A — 1 nuclear system are embodied in the equations (|3.31|) and ( |3.30| ). We solve these 



equations below using a mixed momentum-coordinate space procedure in which the wave 
functions are (p + ,Xj_ | a) = ip a {p + ,Xj_). The values of p + are greater than zero. Thus the 
so-called spectrum condition that positive energy particles have only positive plus-momenta 
is maintained in our mean field approximation. 

An intermediate step is to make an approximation by using coordinate space techniques. 
Here one does not maintain the spectrum condition in an exact manner. Then one can show 
there is a very close relationship (but approximate) between our ip a (x~, x±) and the usual 
solutions to the Dirac equation obtained from the equal-time ET formulation. 

To see this, let's first consider the case where there is no vector potential at all (V^ — ► 0). 
Then multiply Eq ( PUD by 7 + and Eq. fl£3ip by j~. Use 7 ± A T | a) = 7 ± (A+ + A_) | a) = 
7^ | a), and then add the two equations. This gives 

(l°Pa ~ 7 3 (2p + - p-)) Mx~,x ± ) = 2 ( 7 _l ■ p ± + M + g s <f){x-, X± j) V«(aT, x_l). (3.39) 

Convert this to ordinary coordinates using x~ = —2z, so that p + = id + = 2i-J^= — > —ijr- 



' dx dz' 

The operator p + acts as a p 3 operator, and the result (|3.39|) looks like the Dirac equation 



of the equal-time formulation, except for the offending term — p~ multiplying the 7 3 . This 
motivates us to look for a solution of the form ipa(z,x±) = f(z)ip^ T (z,x±), in which f(z) 
is chosen to remove to the offending term. The notation ET refers to the usual equal-time 
solution, because we see that ^ T obeys the usual ET Dirac equation 



provided 



so that 



7° y " 7 ■ P ~ M - g s 0(z, x ± ) J ^ T (z, x ± ) = 0, (3.40) 



f( z ) = e ip « z/2 , (3.41) 

M*,*±) = e ip ^ 2 ^ T (z,x ± ). (3.42) 
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The quantity of interest is ip a (p + , x±) which is expressed as 

i> a (z, x ± )^-= dp + e vp+z MP + , X-0- (3-43) 

V27T J0 

The approximation is that the correct version of ijj a (p + , x±) will have no support for p + < 0, 
but the approximation (|3.43|) does. We can determine this support by examining the inverse 



Fourier transform. This gives 

Mp + ,*±) = -z= / dze-^ + - p ^ 2 >^ T (z,x ± ), (3.44) 

\/ 2tc J— oo 



/2tT 

which is a Fourier transform of the equal-time Dirac wave function at a z-component of 
momentum p + — p^/2. This is not exactly equal to zero when p + is zero or negative, but it 
is very small because p~/2 includes the nucleon mass. The relationship between the term 
p~/ 2 and the binding energy of the level denoted by a is 

p-/2 = M-e a . (3.45) 

Thus the relation (|3.44|) is just the usual equal-time procedure equal-time prescription, 



represented by Eq. (|1 . 7|) , of replacing the kinematic variable p + by the combination of 
dynamical and kinematic variables M — e^ + p 3 for the orbital a: 

p + -> M -4+p 3 . (3.46) 

However, the prescription ( |3.46| ) is dramatically changed when the vector potential is 
included. To see this, multiply Eq. ( |3.31|) by 7" and Eq. (|3.30| ) by 7 + . Then Eq. ( p.39|) 
becomes 

( 7 °(p- - 2g v V°) - 7 3 (2p+ -p~ + 2g v V )) M^,*±) = 

2( lL -p L + M + g s ct ) )i, a {x-,x i _) ) (3.47) 

in which we used V~ = V~ = V°. We again wish to reduce the coefficient of the 7 3 term to 
2p + . This can be done with a new version of the multiplier f(z). We find that the light-front 
wave function is given by 

Mp + ,*±) = ~r= I" dze- i{p+ - p « /2)z e- i9 " A{z ' x ^^ T (z,x ± ), (3.48) 

V 2ll J~oo 



where 



d + A{x~,x ± ) =V°{x-,x ± ), 

POO 

A{z,x±)= dz'V°{z',x±), (3.49) 

J z 



and 



l°(Pa ~ 9vV°)^ T (z, x± ) = ( 1 -p + M + g s 4>)^T{z, x ± ). (3.50) 

The relation ( |3.48| ) tells us that the influence of the vector potential is to remove plus- 
momentum from the nucleons. This removal and enhancement of the nuclear vector meson 
content is the most dramatic result we have. 

How accurate is Eq. ( |3.48|) ? This can only be addressed by solving the problem in a 



manner which respects the spectrum condition. The results show an astonishing agreement 
between the eigenvalues of Eq. ( |3.31|) and those of the equal-time Dirac equation. Thus it 
should be safe to use Eq. ( |3.48| ) for qualitative purposes. 



IV. TECHNICAL ASPECTS 



The solution of the nucleoli and meson field equations are discussed. The reduction of 
Eqs. ( |3.31| ) and (|3.30| ) to a two-dimensional matrix equation is presented here. The new 
numerical technique involving splines is elaborated in App. [A|. 

The nucleon mode equation resulting from the minimization of \{P + + P~) is given 
by the coupled set of equations (|3.31| ) and (|3.30| ). The meson fields <fi and V ± obey the 



equations 

occ 

(-(<9 + ) 2 - d 2 ± + m 2 ) <f>{x-,x±) = -g s J2 $a(x~, x±)i/) a (x-,x±), (4.1) 

a 
occ 

(-(<9+) 2 - d\ + ml) V ± (x-,x ± ) =g v J2 Mx~,*±)^Mx~iX±), ( 4 - 2 ) 



in which ip a (x , x±) = (x , x± \ a). We use the Harindranath- Zhang |13[ representation for 
the Dirac matrices a and (3, which allows us to write Eqs. (|3.31|) and (|3.30|) in 2-component 



form. This representation can be obtained from the standard representation 



by the unitary transformation 



*-*UT- (") 



Hence ip — > Uip and 9 — > U9U* 1 , where 9 is a Dirac matrix in the standard representation. 
In our representation, the matrices of interest are: 

The 4-component wavefunction \ip a ) may now be written in the form 

<*~-^*">=(£:>il«!)- (4 - 6) 

in terms of the 2-component wavefunctions |?/>+) an d iV'a)- Thus the 2-component form of 
Eqs. (|3T31| ) and ( p30|) is 



(p~ - 2g v V~ - 0+)|V+> = (cr ± • (p ± - <?„0 ± A) + tr 3 (M + <? s 0)) |C>, (4.7a) 

za+|^-) = (<r ± ■ (p± - (7„a ± A) + a 3 (M + ^ s 0)) |^+). (4.7b) 

The scalar and vector densities are defined as 
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P S =J2 Ma = E (ft a ^a + ft^t) , ( 4 - 8 ) 

a a 

occ occ 

p± = £ &T^a = E 2(^) f € • (4-9) 

In the nuclear rest frame, p + = p~ = p°, where p° is the usual nucleon density. Hence 
Vj_ = and V~ = V + in this frame. 

A. Angular momentum 

We can write 

(T± ■ p± = CT(+)P(-) + o"(-)P(+), ( 4 -!0) 

where <7(±) = i(cri ± icr 2 ) and 

P(±) =p 1 ±ift = -.- e ± ^(| : ±^. (4.11) 

Here r = |a?j_| and is the azimuthal angle, using cylindrical coordinates. For the nuclear 
physics problems of interest, we anticipate that there is an axis of azimuthal symmetry. 
Hence we can expand the 2-component wavefunctions in eigenstates of angular momentum 
J z , with eigenvalue j z : 

(x- } x±\ij£) = i(x~,r\u^)e i{jz -2^ Xl + (x',r\l^)e i{jz+ 2^ x i 



2 

iu±(x-,r)e iijz ~2 ) ' t ' 
l±(x-,r)e i(jz+ 2^ ' 



(4.12) 



where xi an d X I are the 2-component Pauli spinors. 

2 ~2 

The equations to be solved are then 

(p- - 2g v V+ - t d+) u+ = - (J- + *±I - ig v ^j I- + M*u-, (4.13a) 

p- - 2g v V+ - rd+) lt={^- r - ^ - iQv^j u- - M%, (4.13b) 

id+u- = - (|- + *±1 - ig v ^j £ + M*u+ a , (4.13c) 

id>l- a = (| - ^ - i9«f ) < - M*£. (4.13d) 

The wavefunctions u^ and £, the nucleon effective mass M* = M+g s (f), the vector potential 
V + , and A are all functions of both x~ and r. Eqs. Q4.13 ) have a manifest spin degeneracy 
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under j z — > —j z . Solutions with the same eigenvalue p a are obtained with the corresponding 
replacement 



(4.14) 



Combined with isospin symmetry, we therefore have a manifest fourfold degeneracy of each 
single particle state. The numerical solution to Eqs. (|4.13|) is discussed in App. |A|. 





V. NUCLEAR BINDING ENERGIES 



If these solution to Eqs. (|3.31| ) and (|3.30|) are to have any relevance at all, they should 
respect rotational invariance. The success in achieving this is examined in Tables I and 
II, which give our results for the spectra of 16 and 40 Ca, respectively. Scalar and vector 
meson parameters are taken from Horowitz and Serot [ 24| , and we have assumed isospin 
symmetry. We see that the J z = ±1/2 spectrum contains the eigenvalues of all states, since 
all states must have a J z = ±1/2 component. Furthermore, the essential feature that the 
expected degeneracies among states with different values of J z are reproduced numerically. 
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TABLES 

TABLE I. Comparison of the single particle spectra of 16 in the equal-time (ET) formalism 
(E a - M) with the light-front (LF) method (p~/2 - M). Units are in MeV. 



State a 
Osi/2 
OP3/2 
Opi/2 



ET 



E a -M 
-41.73 
-20.77 
-12.49 



Jz = ±1/2 
-41.73 
-20.79 
-12.51 



LF 



J z = ±3/2 



-20.77 



TABLE II. Comparison of the ET and LF single particle spectra of Ca. 



ET 



LF 



State a 


E a -M 


J z =-- ±1/2 
-55.39 


J z =-- ±3/2 


OSi/2 


-55.40 




0p 3 /2 


-38.90 


-38.91 


-38.90 


OPI/2 


-33.18 


-33.18 




0d 5/2 


-22.75 


-22.76 


-22.75 


lSl/2 


-14.39 


-14.36 




Od.s/2 


-13.87 


-13.88 


-13.89 



J z = ±5/2 



-22.74 



TABLE III. Total plus-momentum per nucleon for 16 0, 40 Ca, 80 Zr, and nuclear matter (NM) 
in MeV. No Coulomb interaction is included here. 



Nucleus 



P^/A 



Pt/A 



P+/A 



P+/A 



16Q 

40 Ca 
so Zr 

NM 



704.7 
672.6 
655.2 
569.0 



6.4 


221.8 


4.7 


253.3 


3.6 


270.2 


0.0 


354.2 



932.9 
930.6 
929.0 
923.2 



22 



The results shown in Tables I — III are obtained using a basis of 20 splines, a box size of 
2L = 24 fm, and 24 Fourier components in the expansion of the wavefunction (see App. |A|). 
This value of L is large enough so that our results do not depend on it, and the number of 
terms in the expression for the density is enough to ensure that the densities are spherically 
symmetric. Another feature is that the spectrum with p + > has no negative energy states, 
so that in using the LF method one is working in a basis of positive energy states only 

The values of p~/2 given in Tables I and II are essentially the same as the single particle 
energies E a of the ET formalism, to within the expected numerical accuracy of our program. 
This equality is not mandated by spherical symmetry alone because the solutions in the 
equal-time framework have non- vanishing components with negative values of p + . 

Table III gives the contributions to the total P + momentum from the nucleons, scalar 
mesons, and vector mesons for 16 0, 40 Ca, and 80 Zr, as well as the nuclear matter limit. 
In the next section we examine in detail the momentum distributions giving rise to these 
expectation values. 

VI. PLUS-MOMENTUM DISTRIBUTIONS AND LEPTON-NUCLEUS DEEP 

INELASTIC SCATTERING 

We discuss the probability that a nucleon, or meson has a momentum p + . In the light- 
front formulation, these distribution functions are determined by the absolute square of the 
ground state wave function. Each distribution is discussed in turn. 



A. Nucleon plus-momentum distribution 

The light-front formulation is very useful for obtaining this observable. The probability 
that we want, Jn{p + )i follows from Eq. ( |4.9| ) as 

occ „ 

/*(?+) = 2 £ d 2 x ± |<p + ,*_l|i£)| 2 , (6.1) 

a J 

with 



oo 



A= dp + f N (p + ), (6.2) 

Jo 

P+= / dp + p + f N (p+). (6.3) 

Jo 

The next step is to define a dimensionless variable y: 

y=p + — = ^, (6.4) 

y y M A Ma 

and a dimensionless distribution fiyfy): 

MV> = ^. (-) 
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The result is shown in Fig. |] for 16 0, 40 Ca, and 80 Zr. The peaks of the distributions range 
from y m 0.72 to y ~ 0.80, whereas the average values (y) are somewhat lower (see Table 
III). The distribution is not symmetric about its average value, as it would be if a simple 
Fermi gas model were used. Both of these effects are caused by the presence of nuclear 
mesons, which carry the remainder of the plus-momentum. 
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FIGURES 
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FIG. 1. Nucleon plus-momentum distribution function, /jv(y), f° r 16 0, 40 Ca, and 80 Zr. Here 
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FIG. 2. Vector meson plus-momentum distribution yf v (y)- In the nuclear matter limit, yfv{y) 
becomes a delta function. 
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B. Scalar meson distribution 

The probability we want is given by 

f s (k+) = fd 2 k ± (W | a\k)a{k) | *). (6.6) 

Using Eqs. ( |3.34| ) and ( |3.10| ) this becomes 



This result which is of the same form as in Ref . [2^| . A final step is to define a dimensionless 
distribution f s (y): 

f,(v) - ^. (6.8) 

The scalar mesons are found to carry less that 1% of the plus-momentum of the nucleus 
(Table III), which is negligible. 

C. Vector meson distribution 

The probability we want is given by 

f v (k + ) = fd 2 k ± £ (* | a^(k,u)a(k,u) | *). (6.9) 

w=l,3 

Using Eqs. ( |3.16| ) and the mean field approximation ( |3.5|) this becomes 

Uk + ) = f ^ 77^-^1 E I Ak, w) | 2 , (6.10) 

J (2vrr (fc 2 + m2) w=13 

in which 

J(fc,w) = f d 3 xe ik - x (ty | $(x)yil)(x) ■ e(Jfe,w) | *). (6.11) 

Using Eq. ( |2.18|) and that only the jj, = u = — term enters, leads to the result that 



d 2 k ± 2k + k\ + ml 
(2vr) 3 (k 2 + mlf k+' 



U k+ ) = / 7^ 77-o 772 X ,L ~ I ^(fe) ! 2 > (6-12) 



another result which is of the same form as in Ref. ||22|| . A final step is to define a dimen- 
sionless distribution f v (y) 

fM = ^j-P. (6-13) 
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The vector mesons carry approximately 30% of the nuclear plus-momentum. The techni- 
cal reason for the difference with the scalar mesons is that the evaluation of a*(k, u)a(k, u) 
counts vector mesons "in the air" and the resulting expression contains polarization vectors 
that give a factor of Ar in Eq. (|6.12| ) which enhances the distribution of vector mesons 
of low k + . The results for the vector meson distribution are shown in Fig. |2|. Clearly as 
the size of the nucleus increases the enhancement of the distribution at lower values of k + 
becomes more evident. In the case of nuclear matter the distribution k + f v (k + ) becomes a 
delta function. 

D. Lepton- nucleus deep inelastic scattering 

It is worthwhile to see how the present results are related to lepton-nucleus deep in- 
elastic scattering experiments. We find that the nucleons carry only about 70% of the 
plus-momentum. The use of our /jv in standard convolution formulae lead to a reduction in 
the nuclear structure function that is far too large (~95% is needed |§]) to account for the 
reduction observed ]l| in the vicinity of x ~ 0.5. The reason for this is that the quantity 
M + g s (f) acts as a nucleon effective mass of about 670 MeV, which is very small. A similar 



difficulty occurs in the (e, e') reaction []25| when the mean field theory is used for the initial 
and final states. The use of a small effective mass and a large vector potential enables a sim- 
ple reproduction of the nuclear spin orbit force flT3,El]. However, effects beyond the mean 



field may lead to a significant effective tensor coupling of the isoscalar vector meson [26 



and to an increased value of the effective mass. Such effects are incorporated in Bruckner 



theory, and a light-front version |27j could be applied to finite nuclei with better success in 
reproducing the data. 

VII. SUMMARY AND DISCUSSION 

The previous Sections present a derivation of a light-front version of mean field theory. 
The necessary technique is to minimize expectation value of the sum P + P + . This leads 
to a new set of coupled equations ( |3.31| ) and ( |3.30|) for the single nucleon modes. These 



depend on the meson fields of Eqs. ( |3.13| ) and ( |3.21| ). 



The most qualitatively startling feature emerging from the derivation is that the meson 
field equations ( |3.13p and ( |3.21| ) are the same as that of the usual theory, except that z of 



the equal-time theory translates to —x~/2 of the light-front version. This can be understood 
in a simple manner by noting that light-front quantization occurs at x + = 0. If one then 
sets z = —t, then x~ = t — z = —1z. However, this simple argument is not really justified, 
because using x ± precludes the use of z and t. A general argument, using the feature that 
a static source in the usual coordinates corresponds to a source moving with a constant 



velocity in light front coordinates, will be presented in a separate paper ||15| . That paper 
also contains a number of solutions of toy models. 

Even though the meson field equations of the light-front and equal-time theories are 
the same, there are substantial and significant differences between the two theories. In our 
treatment, the mesonic fields are treated as quantum field operators. The mean field approx- 
imation is developed by replacing these operators by their expectation values in the complete 
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ground state nuclear wave function. This means that the ground state wave function contains 
Fock terms with mesonic degrees of freedom. We can therefore compute expectation values 
other than that of the field. In particular, we are able to obtain the mesonic momentum 
distributions (Sec. VI). This feature has been absent in standard approaches. 

We obtain an approximate solution ( |3.48| ) of our nucleon mode equation. Our nucleon 



mode functions are approximately a phase factor times the usual equal-time mode functions 
(evaluated at x~ = —2z). This shows that the energy eigenvalues of the two theories should 
have very similar values. But the wave functions are different- the presence of the phase 
factor explicitly shows that the nucleons give up substantial amounts of plus momentum to 
the vector mesons. 

A new numerical technique, discussed in Sec. IV and App. [A|, is introduced to solve 
the coupled nucleon and meson field equations. Our results display the expected 2j a + 1 
degeneracy of the single nucleons levels, and the resulting binding energies are essentially 
the same as for the usual equal-time formulation. This indicates that the approximation 
( ggg ) is valid. 

As discussed in Sec. VI-D, the present results related to lepton-nucleus deep inelastic 
scattering experiments and (e, e') reactions are not consistent with experimental findings. 
This is because, in 40 Ca for example, the nucleons carry only 72% of the plus momentum. 
This is a result of the quantity M + g s <f), which acts as a nucleon effective mass, is very 
small, about 670 MeV. The use of a small effective mass and a large vector potential enables 
a simple reproduction of the nuclear spin orbit force [fP2"p4f . However, effects beyond the 
mean field may lead to a significant effective tensor coupling of the isoscalar vector meson 
|26| and to an increased value of the effective mass. Such effects are incorporated in Bruckner 



theory [p_7 which, for infinite nuclear matter, results in nucleons having about 80-85% of 



the nuclear plus-momentum. A light-front version |27| should be applied to finite nuclei 
with better success in reproducing the data. Another approach could be to use different 
Lagrangians, with non-linear couplings between scalar mesons and the nucleons |]12j], or ones 



in which the coupling is of derivative form [p_HJ : ip^ipducj)- These models are known to have 



significantly smaller magnitudes of the scalar and vector potentials. In particular, in nuclear 
matter vector mesons carry only about 10-15% of the nuclear-plus momentum. Another 
interesting possibility would be to obtain a light-front version of the quark-meson coupling 
model p9| , in which confined quarks interact by exchanging mesons with quarks in other 



nucleons. This model, also has smaller magnitudes of the scalar and vector potentials. 

In any case, these kinds of nuclear physics calculations can be done in a manner in which 
modern nuclear dynamics is respected, boost invariance in the z-direction is preserved, and 
in which the rotational invariance so necessary to understanding the basic features of nuclei 
is maintained. 
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APPENDIX A: NUMERICAL TECHNIQUES 



There are many possible numerical approaches to solving Eqs. Q4.1[j| ). We choose a 
method that is robust, and emphasizes the physical content of the wavef unctions, at the 
expense of being computationally intensive. We begin by making a Fourier expansion of the 
wavefunctions in the variable x~: 



ut(x ,r) 



1 



V2L2 



IX 



£e-^-/ 2 < n (r), 



(Al) 



with a similar expansion for l^ (x~ , r) . Boundary conditions are imposed by constraining 
the system to be in a "box" of a given length in the variable x~ . In the nuclear rest frame, 
x~ = —2z, and so d + = 2<9_ = —d/dz. Hence for — L < z < +L, we write 



u a (z,r) 



1 



V2L2tt 



E^X'nW. 



(A2) 



with {p+ = nq, n = 1, 2, 3, . . .}, and q = ir/L. 

For bound states, the functions u^ n {r) and l^ n {r) are real, and so the scalar and vector 
densities take the form 



P s (z,r) = J2 P s m (r) cos mqz, 



m>0 



P ± (z,r) = J2 P^(r) cos mqz, 



m>0 



with m = 0, 1, 2, . . ., and 



PmiT) 



Pm(r) 



2L2nl + 



I occ 

-L "I Um.O ™ <n 



} m,0 a n 
~''a,n\ r J'a,M-mV ) ~ 'a,nl r J'a,n+mV ) 

2 2 



n\ r ) U a,n+m\ r ) ' U a,nV ) U a,n+mV ) 



occ 



2L 2tt 1 + <L 



2-^/2.^/ \ U a,n\ r ) U a,n+m\ r ) ' ''a,n\ r )''a,n+mv) 



J m,0 a n 

The normalization integral for a nucleus with A nucleons is 



/*oo rL 

A = 2n drr dzp + (z,r) 
Jo J-l 

PCX) 

= 2tx2L \ drrptir) 
Jo 

/■OO occ r 

= 2 drrJ2J:[Kn(r)) 2 + (lin(r)f 



(A3) 

(A4) 



(A5) 
(A6) 



(A7) 
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1. Meson fields 

The equations for the meson fields are solved using Green function methods. We illustrate 
this for the vector field V + , with results for <p following by analogy. Starting with 

( d 2 d 2 1 d \ 

+ m 2 v )V + (z,r) = g v p + (z,r), (A8) 



dz 2 dr 2 r dr 



we expand V + (z,r) in the same form as the density p + (z,r), Eq. (A4): 



V + (z,r) = y^ V^ (r) cos mqz. (A9) 



The functions V^(r) satisfy 

/ d 2 1 



\ dr 2 r dr 
and their solution may be written as 



-7T + < + m V V+(r) = g v pt(r), (A10) 



■^ 



V+(r) = g v / dr'r>G(r,r')pi(r>). (All) 

Jo 

The Green function is 

G(r, r') = I (m* v r)K (m* v r')6(r' - r) + I (m* v r)K (m* v r')6(r - r'). (A12) 



We have introduced the definition m* = Jm 2 + m 2 q 2 } and Jo and K are modified cylindrical 
Bessel functions of zeroth order. The meson fields are computed numerically from Eq. ( [All ) 
by an outward and an inward integration. 

2. Solution of nucleon equation 

To streamline the notation, we drop the explicit dependence on the single particle label 
a in this section. Equation ( [4.13[ ) can be rewritten in the form of a 2 x 2 matrix equation 



^^)^^-^'^T=kl H \(^)' (M3) 



with the constrained subsidiary relation 



(z,r\u~)\_ I n nz,r\u+)\ (AM) 



\(z,r\r) J " (-i£)I \(zAi + ) 
Here I is the 2x2 identity matrix, and 
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n 



M* 

D 2 -ig v § 



Di+ig v § 



-M* 



V or r 



A 



'd j z -\ 
dr r 



If we take N Fourier components n =- 1, 2, 3, 



, N in the expansion of Eq. 

u + (z,r) = (z,r\u + ) and l + (z,r) = (z, r\l + ) have the matrix representation 



/ <{r) \ 
ut{r) 



lt{r) 



(A15) 

(A16) 

(A17) 
in z, then 

(A18) 



where u^(r) = (p^,r\u + ) and l^(r) = (p^,r\l + ). 

Equation (|A13|) becomes a 2N x 2N matrix equation. Matrix elements of the N x N 
sub-blocks are determined from the integrals 



V + (r) 



(nn') 



<p+|V+(z,r)b+> 



1 r L 
2LJ-L 

1 + $m,0 



dze i{p n'- p ^ z V + (z,r), 

V^( r ) S \n-n'\,m- 



Similarly, 



[ M *( r )](nn') = MS n,n> + 7T^9sM r )^\n-n% 



'-if 

az 


(nn 1 ) 


= PtKn' 


\P^\(nn') 


— Di5 n<n 


[D 2 


(nn') 


= D20 n ,n 



[iHr)] (nn , 



(nn 1 ) 



■jV+(r) + (-!)"%» 

Pn ~ Pn' 



\n— n'\,m 



m^O. 



(A19) 
(A20) 
(A21) 

(A22) 

(A23) 

(A24) 
(A25) 

(A26) 



The last relation comes from the definition — §j = V + . Using integration by parts, and a 
careful treatment of surface terms, gives matrix elements in the form ( |A26|) . 

The problem has now been reduced to an eigenvalue problem involving 2N coupled 
differential equations in the variable r. To solve this, we make a further expansion of u^(r) 
and ln(r) in a finite basis of B-splines of degree k pO| , |3~T|j : 



M 



U n( r ) =Y, a ti B 



(*), 



N 



lt(r) = Y,PnA k \r). 



(A27) 
(A28) 



t=i 
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The B-splines {B 



(k) 



, A/"} are polynomials of degree k — 1 spanning a domain of 



equally spaced knots {tj, i = 1, . . . , M + k} in r. They are smooth "local" functions that are 
non-zero only on the interval ti < r < ti + k- This basis forms an accurate, but non-orthogonal 
set. Hence the overlap matrix 



S ij= ^ drrB\ k \r)Bf\r) 

•J u 



(A29) 



is non-diagonal. It is diagonally banded, however, since B\ (r) and Bj (r) are functions 
that only overlap if \i — j\ < k — 1. The property of being diagonally banded also applies to 
the matrix elements of other operators. 

We now have a matrix equation with the 2x2 block structure 



2g v V + -i£ ®I + U 



1 



-n 



a 



^ 



P"[l]®Ji p j, 



(A30) 



where <S> denotes an outer product of matrices. a + = {o^} and f3 + = {Pni} are column 
vectors of length (N x A/"), and the (N x A/") x (N x A/") sub-blocks have matrix elements 



V 



V + (r) 



drrBf\r)Bf\r 

(nn'),(ij) Jo J 

[W] ta n„ fl = rdrrBt\r)Bf'(r)\M-(r)] lm , 



'■I 



{nri),(ij) 







(nn 1 ) 
(nn') ' 



Pn On,n' '&ij 



W\(nn>),(ij) ~ °n,n'&ij 



[Di] {Mj) = K« i drrB?\r) [ ~B?>(r) - 



d rW/^ i z + I 



dr 



B^(r) 



[D 

.dA 

dr 



6n, n , I drrB?\r) \^-Bf\r 



2l(nn'),(y) '•"■" / Q "" ' "% v >\dr~ J 



W/s 3z ~ 2 p (fe) 



5f ; (r) 



"^^(rS^WsfCr))^^)]^, 



(nn'),(ij) 

The last relation follows from using integration by parts. 



(A31) 

(A32) 

(A33) 
(A34) 

(A35) 
(A36) 

(A37) 



3. Numerical methods 



For our numerical calculations we use L = 12 fm and N = 24 Fourier components, and 
set the number of terms in the cosine expansion of the densities to be N/2, or 12. We choose 
k = 5 for the degree of the B-splines, Af = 20 B-splines in the expansion in r, and take 
< r < L. The integrals over r are performed using Gaussian integration between knots, 
which gives exact results for the matrix elements Sy, Eq. ( |A29|) . 

The matrix eigenvalue problem Eq. (|A30| ) is of the form Ax = XBx, where A and B 
are real, symmetric matrices. In our problem, A and B are diagonally banded, and there 
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are efficient EISPACK routines that take advantage of this [32]. Cholesky decomposition is 
used to efficiently compute the matrix 



-£ 



-n, 



(A38) 



which is needed both to determine u and I as well as to construct the matrix in Eq. ( |A30| ) . 



APPENDIX B: MOMENTUM DISTRIBUTIONS 



In momentum space the meson field equations become 



V+(k + ,k A 



9s 



k 2 + m 2 
9v 



p s {k + ,k A 



P + (k + ,k, 



k 2 + m 2 

with k 2 = (k + ) 2 + k\, and the convention 

V + (k + ,k ± )= fd 2 x ± J dze- ik± - x± - ik+z V + (z,x x ). 

The scalar meson momentum distribution from Eq. ( |6.7| ) can be rewritten as 

2k + r d 2 k, 



fsik" 



9l 



2vr J (2tt) 2 (k 2 + m 2 ) 2 



\Ps(k + ,k ± )\< 



2k+ r d 2 k 



2tt J (2?r) : 



\<f>(k + , k x )\ : 



Then 



2 -f I d 2 x ± \<P(k\x ± )\ 2 . 



dk + k + f s (k- 





1 r°° dk + 

2 J-oo 27T 

3„ /a+J./~\Y^ 



2(k + ) 2 fd 2 x ± \(j)(k + ,x ± )\' 



d z x(d + (p(x)y 

(T s ++ ). 



(Bl) 
(B2) 

(B3) 

(B4) 
(B5) 
(B6) 

(B7) 
(B8) 

(B9) 
(BIO) 



The vector meson momentum distribution is a little more complicated. Starting from 
Eq. (JU2D , 



Uk + ) 



2k + r d 2 k 



■'•■ hi +m V(^ + ,fc±)i 2 



2tt J (2tt) 2 (k 2 + ml) 2 (k+) 2 
2k + r d 2 k ± g v k 2 + m\ - (k + ) 2 
~2~7 J (2vr) 2 k 2 + m 2 v {k+)~ 2 



V + (k + ,k ± )p + (k+,k_ 



d 2 k< 



2ttJ (2vr) 2 [k+ 



1 



—g v V + (k + , k ± )p + (k + , k ± ) - k + \V + (k + , kj 



(Bll) 

(B12) 
(B13) 
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Then 

Pi 



CLrv rv Ty I rC 



dk 



2 J-oo 2lT 

I** 



2 / d 2 x ± [g v V + {k + ,x±)p + (k + ,x x ) - (k + ) 2 \V + {k + ,x ± ] 

g v V + (x)p + (x) - (d + V + (x)f 



(B14) 

(B15) 

(B16) 

(B17) 



Clearly f v {k + ) is singular at k + = 0, so we plot k + f v {k + ) instead. 

Momentum distributions involve integrals over x±, or equivalently over k±, so one really 
only needs to Fourier transform V + (z, x±) in z. If we define 



V + (k + ,x ± ) 



dze' lk z V + (z,x ± ), 



then for fc+ = mq, q = n/L, it follows from the definition Eq. (|M|) that 

V+(k+,x ± ) = 2L 1 -±^V+(r), 



(B18) 



(B19) 



with a similar result for p + (fc+, x±), 0(fc+,ccj_), and p s (fc+, cc_i_). Hence we calculate the 
momentum distributions from the expressions 



+ + 2(mg) 2 (l + <5 m , ) 2 o /■» j 

k rnfs(k^) = --* — 2vr 2L / drr(j> m (r)) 

q 4 Jo 



^mJv{^n 



2 ( 1 + ^2Lf^ 



^CWCW - {mYK(r)y 



q 4 jo 

The nucleon momentum distribution for p^ = nq is given by 

2 roc 



(B20) 
(B21) 



MfV, 



g Jo 



drrJ2[(u+ n (r)) 2 +(lt n (r)y 



(B22) 



so that A = f °°dp+f N (p+) « g£ n /*(?+), and P+ = / °° dp+p+f N (p+) « <zEnP£/*(#)- 
We interpolate between the discrete values of fc+ and p+ to produce the plots of Figs. 1 and 
2. 
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